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We study the probability that two directed polymers in the same random potential do not inter¬ 
sect. We use the replica method to map the problem onto the attractive Lieb-Liniger model with 
generalized statistics between particles. Employing both the Nested Bethe Ansatz and known for¬ 
mula from Mac Donald processes, we obtain analytical expressions for the first few moments of this 
probability, and compare them to a numerical simulation of a discrete model at high-temperature. 

From these observations, several large time properties of the non-crossing probabilities are conjec¬ 
tured. Extensions of our formalism to more general observables are discussed. 


Introduction. — Recently there was considerable 
progress in calculating the free energy, and its fluctua¬ 
tions, for directed polymers, or directed paths, in random 
media. This problem arises in a variety of fields, includ¬ 
ing: optimization and glasses [T], vortex lines in super¬ 
conductors [2], domain walls in magnets [3], disordered 
conductors 0, Burgers equation in fluid mechanics [5], 
exploration-exploitation tradeoff in population dynamics 
and economics [5] and in biophysics [71 [H]. Moreover, an 
exact mapping connects the DP in 1 + c? dimension to the 
Kardar-Parisi-Zhang (KPZ) equation in dimension d, 
which, in d = 1, is at the center of an amazingly rich 
universality class, including discrete growth and particle 
transport models, with surprising connections in mathe¬ 
matics to random permutations and random matrices. 

Two very different methods led to exact solutions: one 
based on the limit of discrete lattices, e.g. particle models 
such as q-TASEP, often yielding rigorous results [siiini- 
[iiiin]; the other one based on replica, a standard ap¬ 
proach in the physics of disordered systems m. and 
the mapping to a continuum quantum integrable sys¬ 
tem, solvable by Bethe-ansatz [a HU Ha EQ]. The cal¬ 
culation of the n-th moment of the DP partition sum 
is reduced to the time-evolution of a n-particle quantum 
state, determined by the initial conditions. The evolu¬ 
tion is performed with the attractive Lieb-Liniger Hamil¬ 
tonian, whose spectrum is exactly computable |6l l2T] . 
The derivation based on the replica-Bethe-ansatz (RBA) 
involves some guessing and has often anticipated rigor¬ 
ous results from the math community. For instance, for 
the DP with two fixed endpoints, corresponding to the 
droplet initial condition in the KPZ equation, both ap¬ 
proaches obtain the free-energy as a Fredholm determi¬ 
nant, showing convergence at large time to the Tracy- 
Widom distribution for the largest eigenvalue of a ran¬ 
dom matrix lainHis]. 

An outstanding challenge is to extend these methods 
and results to collections of directed paths with hard-core 
repulsion, a difficult problem involving both interaction 
and disorder in a non-perturbative way. It arises in the 
above examples, e.g. populations competition, steps in 
vicinal surfaces or the vortex glass in 2D superconductors 



FIG. 1. (Color online) Paths with fixed endpoints. By ex¬ 
changing the paths after the last intersection, one builds a 
mapping between crossing paths and paths with the final end¬ 
points exchanged. 

|23) . There was progress in that direction in the context 
of vortex arrays [21], within the multilayer PNG growth 
model [21], and the semi-discrete DP hierarchies [ami, 
with emerging connections to the spectrum of random 
matrices. Within the RBA method, in almost all cases up 
to now, only the Id Bose-gas was considered, i.e. initial 
conditions corresponding to a fully symmetric quantum 
state. Here we consider infinite hard-core repulsion, mod¬ 
eled by a non-crossing condition, which requires however 
more general initial conditions. 

The aim of this Letter is to study continuum DP ob¬ 
servables for non-crossing paths. We develop the more 
general nested replica Bethe Ansatz (NRBA), and con¬ 
nect it to another recently developed method [Sj. Here, 
as a first step, we focus on the calculation of crossing 
probabilities, but we expect the potential outcome of the 
method to be broader. 

We introduce the partition function of a directed poly¬ 
mer with fixed endpoints 

Zr,ix-,y\t)= (i) 

J x(0)=x 

in a random potential with white-noise correlations 
r]{x,t)ri{x' ,t') = S{x - x')S{t - f). Then, the probabil¬ 
ity that two polymers with fixed endpoints do not cross 
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in a given realization r] of the potential, is expressed as 


Pr,{xi,x2;yi,y2\t) = 1 - 


Zn{x2;yi\t)Zn{xi,y2\t) 

Zr,{xi]yi\t)Zn{x2;y2\t) 


( 2 ) 


since all paths with at least one intersection can be ob¬ 
tained from paths with yi,y 2 exchanged [T] (see Fig. [^. 
For simplicity, we will consider the random variable de¬ 
fined by the limit of near-coinciding endpoints 

p^{t) = lim ^"'^ = d^dy\nZ^{x;y\t)\j:=o (3) 

£^o 4e^ y=o 

where the last equality, derived from (H), belongs to a 
larger set of relations between non-crossing probabilities 
and the single path free energy m El]. We now present 
a technique to calculate all the moments of Pyit) at ar¬ 
bitrary time t with explicit results for the first few. 

Replica trick and nested-Bethe ansatz. — The aver¬ 
age of products Zn = Zn{xi,yi\t)... ?/„|t) satisfies 

m 


Z„(x;y|t) = {xi...Xn\e .2/„) (4) 

for any integer n, in quantum mechanical notations, 
where bold symbols are shorthand for ordered sets of 
variables and the Lieb-Liniger Hamiltonian reads: 

n 

Hn = -Y,dl.+2c Y, 5{xi-Xj) (5) 

2=1 l<i<j<n 

with c = -c < 0. To use the replica trick we introduce 


0 «.m(t) = lim[( 2 e) ‘^Z^^\e)]^[Zn(0]0\t)]^ 2 m ( 0 ) 

where we set Z^\e) = Z^(e;e|f).Z^(-e;-e|t) - 

Zy{-f,e\t)Zy{e-,-e\t), so that Pn{t)'^ = 0o,m(t). The ad¬ 
vantage of this expression is that for integers n, m with 
n > 2 m, it can be expressed in terms of Q : 

0„,„(t) = lim(2e)-2-(vI/^(e)|e-‘""|4/™(e)) = 

e^O 

V IMP 

where -\-e,e)) .. .Q) and 

a complete set of eigenstates \p) of Hn of energies 
has been inserted with Here Vm is a 

differential operator obtained from the limit e ^ 0 of 
(2e)-'"(4'^(e)|/i), e.g. 2?i = Since 

Hn is integrable by Bethe-ansatz, the eigenstates, with 
eigenvalues = T,]=i take the form 

n 

E ^^Q(x)HQexp[i ( 8 ) 

P,Q<iSn j = l 

where {/ii,...,/i„} is a set of rapidities, <S„ the set of 
n-permutations and ^?q(x) the indicator of the sector 


XQi ^ XQ^... < xq^. However, |4'm(e)) is not a sym¬ 
metric state under the exchange of the coordinates, thus 
the quantum dynamics described by ([^ does not belong 
to the bosonic sector. Nonetheless, it is still possible to 
explicitly determine the eigenstates [57], corresponding 
to different representations of the symmetric group. It is 
enough to choose the vectors Aq, for all fixed permu¬ 
tation P, inside an irreducible representation of Sn- The 
relevant case for us is the representation corresponding to 
a two-rows Young diagram ^ = (n - m, m), where we de¬ 
note a diagram as the decreasing sequence of row lengths 
[^. For instance, for n = 8 and m = 3 we have 


(5,3) 


1 

3 

5 

7|9 

2 

4 

6 



(9) 


and the filling indicates antisymmetric wave- 
functions under the exchange of coordinates 
xi X 2 ,X 3 X 4 ,X 5 xq, which are in the sym¬ 
metry class selected by the action of 'Drn= 3 - These 
representations can be built explicitly as the Hilbert 
space of an integrable spin-1/2 chain with n sites 
restricted to the sector with m down spins. Then the 
eigenstates of Hn on a ring of length L are obtained 
diagonalizing simultaneously the spin-model. This leads 
to the so-called nested-Bethe-ansatz (NBA) equations 


^ Aab - ic ^ A Ag - - zc/2 

b=l Aab j-i Xa ~ f-lj + ic/2 

bta 

P-jk + ic ^ .pj- Pj - Ag - m/2 ^ 

k=l pjk ~ *C n=l Pj ~ + ic/2 

k*j 


(10a) 

(10b) 


where Pap = Pa - Pp and same for the A’s, the auxil¬ 
iary rapidities on the spin chain that impose the appro¬ 
priate symmetry to the wave-function. Solutions of (10) 
provide the eigenstates of ([^ in the appropriate sym¬ 
metry class and the wave-functions are obtained setting 
Aq = ((^q|w(Hm))) with 


n 

Km)))= ^ a(a|M)cT_“Y..a!-|+)) . (11) 

a\ = l 


Here, {Pp)j = pp^ and |...)) indicates states in the 
auxiliary spin space, is the lowering spin operator 
at site a, a = 1 ,..., to , acting on the reference state 
1 + )) = |t...T)). The vector of states I'I'q)) is fixed by 
the filling of ^ and performs the unitary mapping be¬ 
tween the spin-chain representation and a particular rep¬ 
resentation of shape (n- to,to), such that the exchange 
of two-spins is mapped into the exchange of two parti¬ 
cles. Here = ^2°/Hp^ accounts for the bosonic phase 
scattering with e nj<if(pij), ^2° = and 
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f{u) = ul{u-ic), while 

a(a|/x) = synij, n(l+ f] (A;|/j) 

.k<l \ ■^Ik / 1=1 


i(u|m) = 


■ Mb - ^c/2 

- icl2 u- Hb + icl2 


n 


and sym_^[VF(A)] = Efl W{R\)lm\ is the symmetriza- 
tion of W{\) over the variables A. 

Average of p-qit). — We now consider the case 
m = 1 which selects the subspace of wave-functions 
= X 2 ) = 0. Then, the wave-function in Q re¬ 
mains continuous, even after the action of Pi, and 
we can average over different orderings of the coordi¬ 
nates X. Hence Viif^{Q) = ^ Z_p,q di{Q~^Pfj,)AQ, where 
djnifJ-) = (-*)"*Pme*^'’^|x=o- We then obtain [50] 

-^di(Q-V)l'I^Q)) = N/^E(Ma-A)^-l+)> (12) 

Q a=l 


with Z = [(n- l)n!] ^ ensuring normalization and jl = 
'LbiJ'bln. It leads to \'Di'if^{Q)\^ = n(n-2)!|H^p|P(A,/x)p, 
where T{X,n) = sym^ EaC/^a - /t)/Ca(A|/x)] and we 
note Aq=i = A. The sum over all solutions of (10) must 
then be performed according to ([^, a formidable task 
in general. However, (10b) simplifies dramatically when 
L ^ 00 . For c > 0, the n rapidities pi,...,pn, are or¬ 
ganized in Us bound states, each composed by mj > 1 
particles, with E"=i Wj = n. The rapidities inside a 
bound state follow a regular pattern in the complex plane 
= kj + ^{rrij + 1 - 2a) + named string. Here 
a = l,...,TOj labels the rapidities inside the string and 
are exponentially small for large L. A study of (10) 
reveals that, at variance with the bosonic case, not all 
string configurations are actually allowed, consistently 
with the symmetry of the wave-function [9|. For those 
allowed, using [3IS], we obtain their norm |34j as: 


<I>(k, m) 


4c^ 

c”$(k,m) + 4(/x; - A)2 ’ 

„ (kj - kj')^ + c^(mj - TOj')^/4 
i<j<j'<n, (% - + c 2 (toj + mjv) 2/4 


(13) 

(14) 


For each configuration of rapidities following the string 
ansatz, a multiplet of eigenstates is given by the set 
{A(^\ ..., of solutions of (10a), i.e. Q(A) s 

P+(A)/P_(A) = 1, where P±(A) = Y[i{X- fj,i±ic/2). These 
values cannot be determined analytically for general n, 
however the sum over them can be performed using 
residue theorem 


27rj Q{z) - 1 


where w{z) is any analytic function inside the contour C, 
which encircles all the solutions A^*^ and no other singu¬ 
larity of the integrand. Equivalently, the integral can be 


computed taking the poles outside the contour, which in 
the case w{z) = , are given by Zk = Pk - *c/2 

[51] . The sum can then be performed analytically. More¬ 
over, for P ->• 00 , string momenta become free and we can 
replace 'ffkj rnjLj which leads to 


Bn,l(t) = 



E 


nic‘ 




E 


= 1 n.!(27rc)"» 

rik 

-$(k, m)A„_i(k, m) 


(15) 


with A„ i = (n{n- l))~^h 2 . Here, (mi,... indi¬ 

cates sum over all integers m > 1 whose sum equals n 
and we defined hp = ~ P;)^ “ (*c)^(j - l)^. The 

rapidities pj are written as a function of string sizes and 
momenta according to the string ansatz, so that A„_i van¬ 
ishes on the n-strings. The conserved charges of the Lieb- 
Liniger model have been introduced as Ap = Ej=i Pj i ^2 
being the energy. A crucial property of ( |15| ) is that by 
replacing A„_i ->• A„_o = one recovers the formula for 
Zn{t) = Zn{'x. = 0;0|t) = Qn,o{t) as given in [7]. There¬ 
fore, rewriting A„_i in terms of the conserved charges 
and using the statistical tilt symmetry (STS) (see e.g. 
Appendix of m), we obtain 


Qn,i{t) = 


n{n^ - l)c^ 
12 


-®'-s 


Zn{t) . (16) 


This expression is exact for n>2 and allows the analyti¬ 
cal continuation n ->■ 0. In particular, we obtain 

Pr,(t) = lim0„.i(t) = E (17) 

n^O It 

This is in fact the exact result for Pp{t) without dis¬ 
order, i.e. p{x,t) = 0. This remarkable conclusion can 
also be obtained by averaging ([^ and recalling that the 
dependence of the average free energy of a path with re- 
spect to its e ndpoints is entirely fixed by the STS, namely 
\nZp{x\y\t) = h{f) - {x- y)^/(4t), where 


/i(t) =lnZ^(0;0|f) (18) 

is our averaged free energy (and average KPZ height). 

Alternative derivation. — A different approach was 
recently proposed in [8] (remark 5.25) where non¬ 
intersecting paths were also studied. There, it was pro¬ 
posed a multicontour-integral formula associated to a 
partition of n. We identify the partition with a Young- 
diagram and for the two-row case of our interest, it can 
be put in the form 


Bn,m (t) 


J_ r ^... r 

2 ™ J 27r J 2tt 


^ n fM 

\l<k<j<n 


m 

n d{z2q-l,2q) 

9=1 


(19) 
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where Zkj = Zk- Zj, h{u) = u{u - ic) and the integration 
contours are parallel to the real axis with an imaginary 
part Cj for Zj satisfying Cj+i > Cj + c. Shifting back all 
the contours to the real axis, we encounter many poles 
whose residues reduce to integrals with a smaller num¬ 
ber of integration variables. This expansion can then be 
organized to reproduce the one based on strings in (151, 
with replaced by [9] 






n^l h{li2q-l,2q) 
ni<j<fe<n / if^kj ) 


( 20 ) 


and again the given by the string ansatz. Interestingly, 
is always a polynomial in the /i’s of degree 2m as 
can be seen considering the residue at coinciding points. 


Moreover, Eq. (20) agrees with the result obtained from 


the NBA for m = 1, which gives a completely indepen¬ 


dent check to the proposition in (19). For to > 1, the cal 


culation from NBA becomes more involved but we will 
continue by assuming that ( |2^ retains its validity. 
Higher moments ofprjit). — We now focus on to = 2. 


Upon symmetrization in (20) one obtains 


A„ 


= — 


(n - l)/i 4 - n{n - l)^c^ 


n(n - l)(n - 2)(n - 3) 


( 21 ) 


After tedious calculations, it can be rewritten in terms of 
the conserved charges A„_ 2 ({^p}) [H]- In contrast with 
the TO = 1 case, higher charges, up to p = 4, are in¬ 
volved. It is therefore useful to formally generalize the 
partition function to (t) which is obtained from the 
expression of Zn{t) replacing the imaginary time evolu¬ 
tion with the more general . This is 

the partition function of a Generalized Gibbs ensemble 
(GGE) [35], which we show can be related to a Fred¬ 
holm determinant [3]. Here, we use it as a generating 
function: 0 „, 2 (i) = ^n, 2 ({dp})Z^(t), formally replacing 
Ap ^ dp = and setting /3p -»• 0 at the end. Deriving 
extended STS identities from the invariance pj pj + k 
in Z^(t) , for arbitrary k, we are able to re-express it only 
from the energy A 2 , leading to 

^=-(ld, + ^d^jh(t). ( 22 ) 


Hence the second moment is determined at all times 
from the average free energy h(t) (18). We did not 
find a direct derivation of this remarkable result, and 
it may be a consequence of integrability. At large time 
h(t) -+^(c‘^t)^/^ EEmU] ) so that 


Pp(t)^ - 


12t 


9t5/3 


(23) 


with X 2 = -1.77(1) the mean of the Tracy-Widom GUE 
distribution m- Repeating this procedure for to = 3 we 
can again use Z^{t). Now higher charges are involved 


and the result is expressed as derivatives of a Fredholm 
determinant |9]. It simplifies at large time leading to: 


Pp{ty 




15t 


9t5/3 


(24) 


It is natural to conjecture the leading decay Ppit)"^ 

integer to > 1. However, the knowl¬ 
edge of moments at long-times is not sufficient to recon¬ 
struct the full distribution of p: in view of we further 
surmise that Pri{t) tends to zero (sub) exponentially at 
large t for all but a small fraction ~ l/(c^t) of environ¬ 
ments where typically Pri{t) ~ <?. This is consistent with 
the conjecture |9] Inp^(t) —a(c^t)^'^^ where a = ^-X 2 
is the average gap between the first (X 2 ) and second (X 2 ) 
GUE (scaled) largest eigenvalues, with a « 1.9043 [T5] 
(note that a rs 1.49134 for the hard wall problem [TH]'). 

Comparison with numerics — To check our results, 
we study a discrete directed polymer on a square lattice 
[7] , defined according to the recursion (with integer time 
t running along the diagonal) 

= + (25) 


with Ug £ sampled from the standard normal distribu¬ 
tion. In the high temperature limit /3 « 1, it maps into 
the continuous DP ([^ at c = 1 with x = 4i/3^ and t = 2t/3^ 
[7j. We consider two polymers with initial conditions 
^ti-i ^ ± 1/2 ending at time i dX x = ±1/2. For 

each realization of the Uj f, the non-crossing probability 
p on the lattice is efficiently computed using the image 
method [1]. Comparing with ([^, we deduce p 16p/3'^, 
for /3 0, due to the rescaling of the factor e = 4/3^. 


The numerical results and the analytical predictions (17 
23) are shown in Fig. For the first moment p,,(t) the 


agreement is excellent even at finite temperature, pre¬ 
sumably due to robustness of ([l7| at any time. The nu¬ 


merical check of (23) is more delicate: indeed, the large¬ 
time behavior of the second moment depends strongly 
on temperature and approaches our prediction only for 
/? « 1, see Fig. (right). However, the leading decay 
is found consistent with t~^ down to zero temperature, 
where the polymer paths do not fluctuate thermally and 
for any to > 0: p™ = p. In order to interpolate between 
zero and high temperatures, we conjecture the large time 
behavior of the moments on the lattice: p™ Cm.il3)lt, 
with Cm(/3) - at high temperatures and 

Cm{P °°) = Coo is a constant that we expect to be non- 
universal |3] . This agrees with the intuitive picture of the 
zero temperature deterministic path, weakly perturbed 
by thermal fluctuations. Ghecking the sub-leading terms 
in (23) would require much more intensive numerics. 

Conclusions. — We presented a general formalism to 
calculate the statistics of N mutually avoiding directed 
polymers in a random potential, with explicit results for 
N = 2. Multi-polymer observables are reduced to a com- 


























5 



(b) 


FIG. 2. (Color online) The first (a) and second moment (b) of in the continuum limit are shown vs t for several value 
of /?. Numerical simulations with at least 2 x 10® realizations and time up to £ = 8192. The value of p,, introduced in is 
obtained by Pr/it) = p/(16/3^), with p the probability on the lattice and the time is scaled to t by £ = t/(2/3‘*). Dashed red line: 
analytical predictions, (171 and (231. In (b), the small t expansion p^ is shown with a dashed blue line. Inset in (b): plot 

of p^ £//3'* - P~‘^C 2 {P) vs /? and extrapolation at small /3 = 0 from the best fit with a quadratic function in /3 at fixed t = 10®. It 
shows a finite limit consistent with our prediction 2 ; 32/3 = 10.6(7) (blue circle). 


pact form in terms of conserved quantities of the Lieb- 
Liniger Hamiltonian and expressed at all times by deriva¬ 
tives of a Fredholm determinant, i.e. the GGE partition 
function. As a simplest example we obtained the lowest 
moments of non-crossing probabilities, with an exact re¬ 
lation between the variance and the free-energy, a non 
trivial scaling for the sub-leading part and a pre¬ 
diction of leading behavior for all moments. The full dis¬ 
tribution of the non-crossing probability is under current 
investigation. Going beyond the infinite hard-core repul¬ 
sion remains for the moment elusive, but we are confident 
that further developments of the present method and full 
exploitation of its integrable structure, will allow further 
progress in the elusive interplay between disorder and 
interactions. 
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Supplementary Material for EPAPS 
The crossing probability for directed polymers in random media. 

Here we give additional details about 

• crossing probabilities of paths; 

• the construction of the representations of the symmetric group; 

• the explicit derivation of the m = 1 from the nested-Bethe ansatz; 

• the generalization to arbitrary initial conditions from the expansion of the multi-contour integral formula; 

• the derivation of a Fredholm determinant formula for the generalized generating function. 

• the extension of the STS symmetry and the calculation of the m = 1,2, 3 moments; 

We will conventionally use ^ (and x) fo indicate Young diagrams of size n, with (^*) the size of rows (columns) from 
top to bottom (left to right). We write each diagram as the sequence of decreasing row sizes ^ s (^ 1^2 ■..), Ci ^ ^ 

Equivalently, one can write it as a sequence of decreasing column sizes. There is a correspondence (bijective 
in each case) between these two sets of sequences, and the partitions of the integer n, as ia = 

with when a< a' and the same for . A standard tableau ^ is obtained from the diagram ^ by filling the 

boxes with the integers l,...,n, with the condition that they grow from left/right and top/bottom along each row 
and column. 


CROSSING PROBABILITIES AND THE FREE ENERGY 


Here we show relations, valid in any given disorder realization, between non-crossing probabilities of several paths 
and the free energy of a single path. Consider N directed paths Xi{T) in the continuum, with endpoints 3:^(0) = yi, 
Xi{t) = yi and Xi< X2< --Xm, yi < y 2 < ■■Un- We know from the Karlin-McGregor formula and generalizations [T] that 
the probability that these paths do not cross can be expressed as a determinant: 


IN) 


N 


(xi,...,Xiv;2/i,...,j/jv|t) = n 


1 


2=1 


? 2/2 lb 


det[Zn{xi-yj\t)]NxN 


(SI) 


Consider first N = 2 which is the main focus of this work. Bringing two endpoints close together, i.e. yi ^2 = y^ e/2 and 
dividing by e, equivalently applying the operator 02 {y) = 1(^2 “ Ai)li/ 2 =yi=yi or equivalently for X 12 we obtain: 

P^^\xi,X2;y\t) = dy{\nZrf{x2;y\t) -\nZrj{xi;y\t)) (S2) 

Py^Hx;yi,y 2 \t) = dx{\nZy{x;y 2 \t)-\nZy{x;yi\t)) (S3) 

Bringing points together on both ends we obtain the non-crossing probability for 2 paths both from y to x: 

P^rfHx;y\t) := 02{x)02{y)p^^Hxi,X2;yi,y2\t) = dxdy\nZv{x;y\t) (S4) 

In particular this leads to Eq. ^ for the quantity defined in the text Py{t) := 0|t). Note that the STS symmetry 

also implies 

P^y\xi,x2]y\t) = , P^y\x]yi,y2\t) = ’ P^v\x]y\t) = (S5) 

in addition to the result mentioned in the text. Note that these STS results are exact also for a model with a more 
general noise, provided it has the STS symmetry, e.g. r]{x,t)r]{x',t') = 6{t- t')R{x - x'). 

Similar relations, although more involved, exist for more than 2 paths, N > 2. For instance the non-crossing 
probabilities with one endpoint coinciding, defined as Py^\x;yi, ..yN^t) = ON(x)py^\xi, ..XN]yi, ..y]\[,t) where now 
0{x) = A ni<i<j<Ar(c/! “ di)\x^=x^ admits a simple expression as a determinant: 

(TV)/ |,^ , .vdl~^Zy{x,yj) 

P/ ’{x]yi,..yN\t) = det [——- ^\nxN 

'^ri\X, yj ) 


(S6) 
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the first row of the matrix being the vector (1,1). Such derivatives can then be re-expressed in terms of derivatives 
of the free energy. Taking the second endpoint coinciding by applying Ojv(y) leads to, e.g. for iV = 3: 

pl^\x]y\t) = {x;y\t)da:dypl^^\x;y\t) - {x;y\t)dypl^'> {x-,y\t) + (S7) 

where we recall p^\x;y\t) = d^dyln Zy(x- y\t), and to more complicated relations for higher N. 


Dependence in elastic coefficient and a conjecture for discrete model: 

If, in the continuum model described by Eq. the elasticity term /g is replaced by /g dT^{^Y lii Ed- 


( s'! 

(1) to study the effect of elasticity k, the above result is trivially changed into py {x-,y\t) = 

Let us now consider a discrete DP model on a lattice, e.g. as the one defined in the text. This model does not 
satisfy exact STS anymore. However it can usually be described by an effective elastic constant Keff{T) defined 
from the curvature around a minimum of the average free energy with respect to one end-point position. It may, 
in general depend on the temperature (defined for the lattice model). Hence we can conjecture that the large time 
limit of the observable p defined in the text will be p ~ ^We/HT) discrete model defined in the text 

oo) = 1). 


REPRESENTATIONS OF THE SYMMETRIC GROUP 


The eigenfunctions of in Eq.([^ are contained in £^([0,L]"). Because of the integrability of the model, the 
eigenfunctions can be written as linear combination of plane-waves in each sector xq^ < ... < xq^ as in ([^. Since 
Hn is symmetric under the exchange of coordinates, they can be classified according to the representations of the 
symmetric group. For a fixed P, the vector of components Aq belongs to a vector space V of dimension nl. In this 
space it is naturally defined a representation of the symmetric group called the regular representation, where the 
action of a permutation T 6 is simply the left multiplication and the matrix representation is C{T)q^qi = 5t~^q,Q'. 
In a similar way, it is defined the dual representation based on the right multiplication 'R,{T)q qi = associated 

with the exchange of two coordinates (see below). Note that C{Ti)C{T 2 ) = C{TiT 2 ) while TZ{Ti)Ti{T 2 ) = TZ{T 2 Ti)- 
moreover C{Ti)TZ{T 2 ) = TZ{T 2 )C{Ti). 

Moreover, the requirements for Q to be an eigenstates imposes that [5] 


ACTiP 


- E 

Q'iSr. 


[Eplp.^jQ.Q'^Q' 


(S8) 


where <Ji is the transposition permutation exchanging -f 1 and = f{pba)C.{cfi) + ifipba) - 1)^(1) with 1 
the identical permutation. This ensures that all the vectors Aq for different P can be chosen inside an irreducible 
representation of iS„, which are in one-to-one correspondence with Young diagrams [3]. It is useful to summarize the 
construction. Given a standard tableau we associate two subgroups of Sn- -d(^) and B{f) defined respectively as 
the two subgroups that preserve respectively rows and columns, i.e. .A(^) are all permutations within rows, and B{f) 
within columns. Two elements of the group algebra of Sn are then defined as 

aiO= E H0= E (-1)""^ (S9) 

PiAU) P6B(U 

from which we define the Young symmetrizer 

ciO = aimo = E cpP . (SIO) 

PcS„ 


The coefficients cp are integer numbers obtained from ( |S9[ ). The Young symmetrizer projects onto = TZ{c{^))V, 
an irreducible representation under the left-action £(iS„)V^, where TZ has been extended on the group algebra by 
linearity e.g. wave-functions built using the vectors Aq in are anti-symmetric on 

the variables inside each column. Indeed for any permutation T € S^. 

V’m(Ex)= Y, exp[iXla;Q.MP.] (SH) 

P,Q^Sn 'i- 

using that {QTK)j = {Tx.)q^ = xtq. and performing the relabeling Q QT~^. Now, for T e B{fj) one finds that: 

V'm('Ex) = (-1)'"'^ Y ^^Q(x)^Qexp[l^a:Q,Aip.] (SI2) 

P,Q^Sn i 
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In the last equality, we used that for any vector v e V^, v = TZ{b{^))v' for some v' e V, and therefore TZ{T)v = 

nnoTw = = {-ir-Tzm))v' = i-iy-v. 

An alternative path is to build the irreducible representations of using the Hilbert space G of n spin 1/2. The 
action G{T) of a permutation T is defined as the permutation of the spins 


G{T) |ll . . . in} = |*Ti ■ • ■ *T„)) 


(S13) 


where 4 € the two eigenstates of the k-th spin along the z-direction. The sub-space of highest weights, 

i.e. annihilated by = Efcs/;, with fixed total magnetization = (n - 2m)/2, defines the irreducible 

representation corresponding to the two-rows diagram ^ = (n- to,to). This procedure can be extended to the general 
case and allows deriving the equations for the rapidities fi, together with a hierarchy of auxiliary rapidities for each 
row of called nested-Bethe-Ansatz equations [4]. For the two-rows diagrams, we get (10) in the text. Although these 
equations only depend on the diagram wave-functions depend on the tableau For instance, for n = 3 and to = 1, 
we have two possible tableaux 




(S14) 


All the possible standard tableaux correspond to different copies of a given irreducible representation inside the 
regular one and therefore to different multiplets of wave-functions. For a general two-row diagram the dimension of 
the irreducible representation is 




= d 





n - 2m + 1 
n - TO + 1 


) 


(S15) 


and its multiplicity inside the regular representation is also equal to d^ [3]. The spin-representation can be mapped 
into each of these copies by a linear mapping dr; H ^ G, with components 'b(r') = Y,q vq l^tg)), for any vector v eV 
of components vq. A complete characterization of this mapping is obtained requiring that, for each choice ^ for the 
same 

• 'b(V'j) s G^, unitarily so that, for any w 

(S16) 

Q 


• for any w e 1^: 


'i>i£{T)w) = G{T)'i>{w) 


(S17) 


• we fix the action of 4' on the space orthogonal to by Y,q I^'q)) ((^q| = with Hj the projector of G onto 
With these definitions, we take 

A^ = A^^A{'^QHPt^)} ■ (S18) 

In practice, the explicit form of dt is not needed, since the differential operator 2?^ gives a non-vanishing result only 
for one Young tableau corresponding to the filling described in (|^. In order to show this, we now give the explicit 
derivation of Eq.(12) in the main text. The vector of components Vq^ = dni{Q~^fi) belongs to and we only need 
the image |to)) s = T.q dm(Q~^fJ-) I'I'q)) of this vector in the spin representation G^^. To determine it, we 

expand it as 

n 

\m))= £ qai,....a„(H)cr!Y..Cr!™| + )) (S19) 

ai = l 


where (7ai,...,a„( m) is a symmetric tensor obtained as a linear combination of the Its components are therefore 

homogeneous polynomials in the of degree TOBecause of (S17), under the mapping 4', a permutation of the rapidities 
/X corresponds to a permutation of spins, we deduce the expansion 


9ai,...,a„ (a*) - ^ 9fc(A*)9m-fc(Mai i • ■ • i ) 

k=0 


(S20) 
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where are symmetric polynomials of degree k. Moreover, from \m)) = 0, we have 


Qaa 2 ...am 0 ■ 


a=l 

a^a2,....,<in 


These conditions are sufficient to determine the state |m)) and in particular for m = 1, they lead to 

qa = \^{Ha - A) ! 


(S21) 


(S22) 


where A = ^ Ea=i Ma and, since 'I' satisfies (|S16|), Z can be fixed equating = ((nr|m)), recovering 


NESTED BETHE ANSATZ APPROACH 


Norm of a state: general form 

The expression of the norm for a given set of and A solutions of ( [To| on a circle of length L takes the form 


f lV'M(a:)pdx= ^ (E^q(^q)*) /" 

P,P’<iS„ \q I 


(S23) 


Now the sum over Q can be performed employing Since I'I'q)) (('I'qI = IIj, it acts as the identity inside 
and we have 


^ f ((a;(PV)|a;(P/x))) . 

p p, J 


(S24) 


This last expression is the norm of a state of a two-components Lieb-Liniger model. This integral can be computed 
exactly in the framework of algebraic Bethe-ansatz [5], leading to 


||V'^ir = c'"B(A)det( 

\ ^xx 

where the matrices e nx G^\ enx Gxx ^ mx m are given as 


{Gfj,fj,)jk - Sjk{L QjaKjk 

a=l 1=1 


{Gfj,x)ja - {Gxfi)aj - Q_ 




{G,,)at = Sat IE Qja - E ^a' | + 


a'=l 


and we introduced the notation 


B{\) E n 


l<a<b<m \ ^ab I 


(525) 

(526) 

(527) 

(528) 

(529) 


Qja — 


- XaY 




2c 


c2 + 


K 


G) = 

ab 


2c 


+ (Aa - XbY 


(S30) 


String states and their norms 

Examination of the NBA equations, and consistency with the nested contour integral method (see section below), 
indicate that at large L the rapidities of the eigenstates fi are arranged in a set of strings, as is the case for bosons: 


a , 


+ 1 — 2q) + 


(S31) 
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where a = 1,... j = 1,... ,ns and Us is the number of strings, while mj > 1 is the size of the j-th string (called 
an m^-string). For rrij = 1, = kk is real (called a 1-string). The Sj are the deviations from the string hypothesis, in 

general different from the bosonic ones, but that we can still assume to be exponentially small for large L. 

An important difference with the bosonic case is that now some of the strings are missing. Let us call the factor 


containing Ao in (10bI. In the bosonic case m = 0 this factor is set to unity and one recovers the usual BA equation: 
from its poles and zeros one sees that the strings (S31) are the solutions. For m > 1 the factor 11^, upon substitution 


of \a from solving (10a), will sometimes cancel some of these poles and zeros and some strings will be missing. Let 


us show two explicit examples for m = 1, in which case there is a single Aa = A. 


Consider n = 2, to = 1. Then the solution of (10a) is Ai s A = + fJ- 2 ) and (10b) becomes simply = 1. 


The 2-string states are missing, the only states contain two 1-strings, with, in fact free particle momenta quantization. 
This is expected from the antisymmetry of the corresponding Young diagram which make the two particles fermions 
unaffected by the 6 interaction. 

Consider now n = 3, m = 1. One looks, with no loss of generality, for solutions such that + /i 2 + /is = 0, and such 


that /is is real and fi 2 = fJ-l- From considering the modulus of (10b) for j = 3 one finds that A must be real. On the 
other hand there are now two solutions to (10a) A = s- 4(/ii + ^ 2 )^ - 4/ii/i2. At large L one then finds two 


types of solutions of (10): either (i) three 1-strings with all fij real and the usual quasi-free quantization conditions. 
Or (ii) /ii _2 = |(fc ± is), s > 0, and /is = -k. The r.h.s. of (10b) for j = 1 vanishes at large L, hence we look for zeros of 
the l.h.s. There are two such zero es: one for s = c, which leads to the 2-string plus 1-string state, hence allowed here. 
The other zero requires s = 2c and k = 0 which is excluded by the condition that A is real (which implies + 3fc^). 

Hence the 3-string state is missing. 

We will continue by assuming that in all cases this structure remains, and that for allowed string configurations, 
the Aa do not introduce additional singularities in the modified Gaudin matrix displayed above. Then following [0], 
we deduce the following large L limit for 


n n 


1 


• ^ A® 

a=l 




n 


TO,- 


J U=1 


c'"L"=H(A)detGAA 


Now consider the factor 


tnO ) 2 _ rr (Mi h-j) _ -p-r h-j Mt 


and again we insert the rapidities organized in strings as in (S31). To do so we rewrite 


n 


a ,,a' 


This product splits into two contributions. The products inside same string j that we label Oj and the products from 
two different strings j + j' that we label Ojj'. These two contributions can be written as 


0, = n 


M--M- 


a^a r^j 


M? - M? - 




nr/a'=i(*c(a-a')) 


<5“) / (*c(a - a' - 1)) / to, 

ata +1 


(n 


1 


Ojf =n n 


Z kj - kf + 


3 3 ^ ^ ^ a^a '''' ^^ ^ 3 3 

ata +1 

^ .1^, . m,/ . . , 


2 — + ic{a' - a) PI kj - kj> + 2 


— + ic(a' - 1) 


=1 a'=i kj - kj' +-5— — + icya' - a - 1 ) a'=i kj - kj' -^ —— + ic{a' - 1) 


(532) 

(533) 


So that we have 


kf + -ica\( k. - + 


ic{mj-m-r) . / ^ N, \ / 


- ic{a - 1) 


-It 7 icimj-m./) . I 1 7 7 icimj+rrij/) . , . >. 

a=i \ kf - kj' + -^— - tea ) \ kj - kj' + -^— - ic{a - 1) _ 


{kj - kj,y + 


2 




(S34) 


and finally 


Kr = n o, n On'Or, = n—n m). 


\j<j 


. a S^-Sr 


(S35) 
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It follows that the following ratio is finite 

(^0)2 c"-—ci>(k,m)n,^ 

IIV'mIP " i"=5(A)detGAA 


(S36) 


Note that, in the bosonic case m = 0, in order to enforce the condition Eg I'I'q)) (('I'qI = one has to choose 
|vI/q)) = Inserting this expression in (S18), we see that with these conventions This is the origin 

of the missing n\ factor in (S36), compared for instance to [7]. 

It is easy to check the formula (S36) in the case of ris = n 1-strings, i.e. all real. Then in (S24) one finds that at 
large L the term P' = P gives the leading contribution, hence: 


L 


= ^EI^qI 


P.Q 


p 


P | 2 | 




(S37) 


For simplicity let us consider m = 1. Then = Ea=i «^o(A|/x)cr“^ |+)). Since all are real (as well as A, as discussed 

above) we find: 




in agreement with (S361 and the formula given in the text (for rrij = 1 and Ug = n). 


(S38) 


Explicit derivation of m = 1 case: re-summation over the spin-rapidities 


The action of the differential operator T>i on the wave-function in ([^, after averaging on the different sectors 
ccQi < XQj ..., can be written explicitly using © in the text and leads to 




n\ 


V2 


di(Q-iPM) 


The sum over Q can be performed using (12), leading to 


7^V'/x,a(0) = - A)Ka(A,PM) = iVZn\ 


a=l 


Mi - Mj 


iM<n Mi - Mj - ^c 




where takes the explicit form 


^(a,m) = ^E 


Yy MPj - MP, - 

l<j<l<n ^^Pl ~ MPj 


■^ / p\nd=a+i(A MPd fc/2) 

(.MP„ - Ml ^ 

nd=a(A-MPd+ *c/2) 


Therefore, the norm square is equal to 


|2?V'M.A(0)P = n(n-2)!|F!^P|^(A,M)p 


(S39) 


(S40) 


(S41) 


(S42) 


By taking into account the norm of the wave-function and (S36), we see that the ratio remains finite 

in the L ^ oo limit. In this limit the string momenta kj become arbitrary real number but a multiplet of wave-functions 


is obtained in correspondence of the set of {A*}, solutions of (10a). Specializing (S36) to m = 1 we obtain 

*c-ip(A*)P+(A*)|^(A*,M)p 


I- 

A* 




rvris-n T Tis 


5 IN-1){P+{X*)P!_{X*) - PC(A*)P_(A*)) 


(S43) 


An,l(p) 


where we used that for m = 1, we have the equality 
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and we have B(X) = 1 for m = 1. The sum over the solutions of (10a) can then be replaced by 

dz P(Z)P+(Z)|^(Z,/X)|2 






2Trc P4z)P!_{z) - Pi{z)P-{z) 




dz |P(z,m)|^P_(z) 
27rc P-{z)-P+{z) 


where the contour encloses only the roots A* and no other singularity of the integrand. The integral can also be 
computed by considering the poles outside the contour, i.e. the poles of \P{z,^jl)'^, at the zeros of P+{z): Zb = fj,b-icl2. 
We have therefore calculating the residue 


An,l (/^) — 


1 

-psym 

n - 1 


fll - - ic 


n , , 

^j<i d-i d'j 


5(m) 


(S44) 


where we set 




j<l ^Qi ^Qj 


T^id-a- - fl) 


a,a 
b>a 


Y\.d>a(.f^b f^d Y[d'>a'(.l^b f^Qd') 

Ud>a{fJ.b - fJ'd) nd'>a'(Mfe-MQd -*C) 
d^b 


(S45) 


where in taking the complex modulus square we have used the fact that fi and its complex conjugate are identical up 
to a permutation. 

In order to make more explicit the expression for A„ i(/x), we start noticing that g{fJ,) is a polynomial in the /r’s, 
as 5 (/x) has no singularities. Indeed 

• For p,Q, = At /3 + e, for arbitrary a, 13 we can have singular terms of order up to e~^. For a,P > a, we can possibly 
have a term of order but then the terms with b = a and b = (3 cancel each other. Singularity of order e“^ 
come from the first product term in (S451. But this cancels out between the permutations Q and , where 
Gap is the transposition exchanging indexes q;,/3. 

• For fia = - ic + e, we can have a singularity of order e“^, coming from the last product in ( S45| ), when Qd = a 

and b = f3. Therefore for this to be there, we need {Q~^)p > {Q~^)a (otherwise the zero of the first term would 
cancel the singularity) and < a', which shows that the condition is never realized since {Q~^)a ^ o'. 

Now, with a similar procedure it is easy to see that also A„ i(At) has to be a polynomial. Moreover it is symmetric 
in the /i’s, and homogeneous in the g’s and c, and is not changed under the transformation + z for any z. 

Therefore 

A„.i(At) = + cG(i)(/x) + 

n-1 

where is homogeneous of degree k in the a^’s . Now can be computed setting c ^ 0. In (S45l, one 

sees that the limit c ^ 0 imposes a = Qa' = b and therefore 






Moreover G^^^(At) = 0, being the only possible symmetric polynomial, function of differences, of degree 1. Since 
g{fj-) = 0 for fdj = -ijc, we deduce that G^^ = n(n^ - 1)/12. Finally we get 


I 

A* 


|PV^m,a(0)P 


(u-2)!4>(k,m) rij 


'/X. AI 


nTls-n T 


„ 2 


\^<3 


12 




where in the last expression the fj,’s have to be organized in the string ansatz. 


GENERAL INITIAL CONDITIONS: RESIDUE EXPANSION OF CONTOUR-INTEGRAL 

Suppose now we are dealing with an initial condition, whose properties under exchange can be encoded in the 
Young diagram ^ filled with coordinates x: we require the wave-functions to be antisymmetric when two coordinates 
in a column of ^ are exchanged. In [5] the disorder average of products of partition sums of groups of non-crossing 
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Zl,l 

Z2,l 

Zl,2 

t-" 

Z2,2 



•2^0:,1 

Za + 1,1 


zn-m,i 







Za,2 



FIG. SI. Left — Sketch of the top-bottom/left-right filling of variables for a two-rows Young diagram. The ordering of the 
variables is represented by the arrows. Right — Structure of the poles (black disks) at Za,i + 1 with a > 1 when integrating the 
variables zij, j = 1,2,3 with ^ = (321). 


directed paths was considered. Each such product can be associated to a Young diagram ^ and involves ndO groups, 
each group containing paths mutually non-crossing within each group and with coinciding starting points. Based 
on a similar theorem obtained there for the semi-discrete directed polymer, it was conjectured (remark 5.25 there) 

that for the continuum model the global partition function Z^{T) = Z^c{T,Xa) (see notations there) can be 

written as a multiple contour integral 


ncii) 

2dT)= n 


(27r*)?“^° 




dz, 




n 


l<a<P<nc(0 \l=li = l 



(S46) 


where the integration domain for the variable Zq-j 
C l > (72 + 1 > ... > Cfc s- (fc - 1). Here we set 


is chosen along Cq + IM, i.e. parallel to the imaginary axis, with 


nc{€) C 

F{{z})= n n< 

a=l j=l 


,+JCoZc 


and in particular we assume it to be is an entire function of the variables {z}. For 

= 0, T=2C% (T) ^ en,m 


(S47) 


(S48) 


we recover (19) when ^ = (n- m,m), i.e. a two-rows diagram with = ... = ™ = 2, = 1,..., = 1. To 

see this it is enough to perform the change of variables Zaj ->■ izAjc, where H = 1,..., n, with the ordering depicted in 


Fig. SI (Left). Note that contours change as given in the text, where we have also distinguished contours for variables 


with the same index a (which is immaterial since no poles are encountered when bringing them together). 

We will now show that, as for the bosonic case [5] where all ^“ = 1, this expression is equivalent to the spectral 
expansion Q in the text, into string solutions of the Bethe Ansatz equations. We use the notations of and extend 
the method to the case of a general symmetry. Starting from Zi_i and the contour Ci, we recursively move all the 


contours to the leftmost C„p(^). While displacing the contours, all the residues of the first product in (S46) have to be 


collected: they produce a hierarchy of additional terms with fewer integration variables Za^.^i^, which are the ending 
points of a sequence of residues Zai,ii = Za 2 ,i 2 + 1 = + 2 = ... = + fc - 1. Each of these terms is characterized 

by the length Xa = k oi these sequences, that we arrange in another diagram of size n: y = (xiX 2 • ■ filling (j) of 
the diagram y corresponds to assign the variables Za,i to each box thus fixing the sequence of residues. For instance 
for ^ = (42) and y = (321) one possible filling is: 


Zl,2 

23,1 

24,1 

Zi,i 

22.2 


Z2,l 




zi,2 — zsx + 1 — z^x + 2 
^1,1 = ^2,2 + 1 


(S49) 


and in general we indicate as z^^ ^ the variable in the box at row i and column j of the filling (j). In order to have 
a non-vanishing residue, each index a has to grow strictly along each row. This restricts the possible diagrams y to 
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those satisfying x ^ meaning that either ^ = X ot > Xi^ where i is the smallest index where they differ. Then 
Eq. (S46) is expanded as a sum over all the elements of the set of allowed fillings $(x) of the diagrams x ^ C 


1 


^AT)- n 

Q=1 S • x<£ 

X=1“12“ 




1 r r dz,h 

§ * 


P,Xp 

2'Ki 


X Res0 


nn 

la</3 t,j 


a,i ~ ~ 


('‘T'CiO C 

n 

CK = 1 itj 




(S50) 


The multiplicities ak are the number of indexes a such that Xa = k and avoid the overcounting when a diagram x has 
multiple rows of the same length, i.e. multiple strings of the same length. The expression Res,^[G({z})] entails the 
iterative computation of all the residues at z^f. ^ = + 1 for k = 1,... ,ni^(x) and I = 1,..., x/c - 1, obtaining finally 

a function of the rightmost variables z^^ in each row of (f). Now we consider the relabeling 

('^02,1 ? '^<i*2,2 J • ■ ■ 1 ■^<i>2,X2 ) ivx 1 +X 2 1 yxi+X2-i j ■ ■ ■ 12/xi+i) (S52) 

: (S53) 


which we indicate shortly as Za^i 
also label the final variables wj - 


- y<p*{a,i)'i this defines implicitly, for each (j), the mapping (jj* : {a,i) {1,..., n}. We 

for 1 < j < npi{x)- For instance, for the filling (j) considered in (S49[) 


Zl ,2 

2:3,1 

24,1 


ys 

y 2 

yi = wi 

^1.1 

22,2 


ys 

y 4 = W 2 


2:2,1 



ye = UJS 




(S54) 


In this way, we arrive at 


"c(^) 1 

2dT)= n 7X7 E 


=\ aiW.- 


E /■•■/ n 


_ p 

27Ti 


X ReSx 


nn 




a<l3 i,j yit>*{a,i) y<l>*{l3,j) ^ / \ a=l i^j 


ncii) C \ 

n Yliyric^ ,i)-yr{c‘,j)) n({2/0*}) 


(S55) 


Notice that in this expression we replaced Res^ with Res^^, since when written in the variables y’s, the sequence of 
residues to be computed does not depend on the specific filling (j) but only on the diagram x- Now, we observe that 
if the constraint x ^ C is removed and the sum is extended over arbitrary filling (j) of the variables y’s in Xi the result 


does not change. This because the added terms will have a vanishing residue. For instance, if in (S54) one exchanges 
02,1 and 01,1 there is no associated pole in the first factor. In this case the sum over all possible fillings (j) without 
constraints, amounts to a sum over all permutations of the variables y. Hence, the above expression can be replaced 
by a symmetrization over the variables y, and we arrive at 


"c(4) 1 „| f 

2dT)= n ^ E n ^xResJsym^Wo({y})] . 


(S56) 


with Ho({y}) defined as 


^o({?/})= nn 

\a</3 i,j 


*{a,i) y<f>*{y,j) 


*{a,i) y4>*{P,j) 


I 


_ j rr yk yi I j -p-r yk yi 1 
Ui j/fc-y;-!/lU yk-yi 


n j'^({y0*}) ■ 

, «=1 i*3 / 

\ /"c(5) 5“ \ 

n Yliy^iic. ,i) y0*(aj))(y(/)*(Q,_?) y(p*{o',i) i) j-^({y(?!)*}) ■ (S57) 

/ \ OL=l i>j I 
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Here, in the first line, because of the symmetrization the filling (j) oi ^ is arbitrary. In the second line of (S57), and 


everywhere below, we choose the “canonical” one (j) = (j)^ satisfying + 1) - = 1, and increasing with 

a (see (S61) below for an explicit example), in agreement with the ordering in Fig. SI - note: do not confuse the 


filling (j)^ of the diagram used to define Wq in (S57), with the the filling of x in (S54). The only term singular when 
evaluating the residue in this expression is the first one, which is already symmetric over j/’s, and can be expressed as 
a determinant of a nfj(x) ^ matrix [5] 


ReSv 


■p-r Vk Ul 

Vk-yi-^ 


= det 


\w, + x^-Wj) 


n-Rix) 


(S58) 


The remaining term is not singular upon evaluation of these residues and we just need to insert the variables obtaining 
our main result: 


Proposition: — 


"c(5) 1 I ^ r d?/: / 1 

^d'r)= n ^ E , f - f n y^xdetj—--1 Subx[sym VF^({?/})] . (S59) 

a=i C“! ai!a 2 !-"-7 J p=i 2m \wi + Xi - Wj /^ 

Here the operator Sub^^ requires replacing = yi_i + 1 for for i {1,1 + xi, 1 + Xi + X2, ■ • ■ j 1 + Xi + • • ■ + Xn_R(x)-i} 
and considering the result as a function of Wj - “ 1; ■ ■ ■ i^i?(x)- The function H{({y}) contains the 

non-symmetric part of lTo({y}) and is the only part depending on the choice of the diagram 


W^C({2/}) = 


(n 


Vk-yi-^ 

yk - yi 




0 y<l>*f{oi,j)){y(l>*{a,j) yif>Ua,i) 1) I }) ■ 


(S60) 


a=l i>j 


with (j)^ the “canonical” filling associated to the diagram 


Remarkably, all the other terms remain equal to the bosonic case. Finally, for the two-rows diagram ^ = (n- m, m), 
Eq. (15l and (20l are recovered by the change of variables Wj ->■ ikjjc- (mj - l)/2, j = l,...,ns, and in general 
Pa -*■ *AtA/c for A = l,...,n, where each row of x is identified with a string configuration, so that mj = Xj 
n-s = njj(x). The substitutions in (S48) are finally applied. For instance, for ^ = (42), we have the filling 


Cj}^ = 


^1,1 

^2,1 

2:3,1 

2:4,1 / 

2/1 

2/3 

/ 2/5 

2/6 

^1,2 

^2,2 / 



2/2 

2/4 




(S61) 


and we obtain 


n (y<^;(a.i)-2/0Raj))(2/0Raj)-y0‘(a.i)-l) = [ (y2 " yi ) (Z/l " ^2 " 1)] [(^4 " (j/S " 2/4 " 1) ] = 


Hhi,2)h{y3A) 


(S62) 




To avoid confusion, we remark that Eq. (S61) is used to define in (S60). Only after the symmetrization over y in 
(S59) has been performed, one must carry on the replacement Sub^ corresponding to each diagram x, as indicated 
above, which is equivalent to injecting the string eigenstates. 


In general, the determinant in (S59) produces the term $(k, m) in (15), while the first and second factor in (S60) 
give respectively the denominator and the numerator in (20). 


GGE PARTITION FUNCTION 

We derive a Fredholm determinant form for the generalized generating function, defined as 

J^(u,t)=^^^Z^(t). (S63) 

n\ 
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where, in this work, we introduced the GGE partition function as the following sum restricted over bosonic states (i.e. 
fully symmetric) 


z^it)= E 

fl, bosonic 


lEfM!!.e-‘^2+Ep>i/3p^p ^ ^ 


n\c 


177^ 


= 1 ns\{2-KcY^ = 


Us /- 

s n/ 

.1 = 1 


+ CX3 ^^2 + Ep>l 


$(k, m). (S64) 


the charges being defined as Ap = Y,j Mp the second equality uses the decomposition over string states valid in the 
large L limit. It is thus a generalization of the bosonic partition sum of the fixed endpoint directed polymer studied 
in [71 . Her e it a ppears naturally in the non-symmetric problem, and its calculation generalizes the one of Ref. |7] . 

Note that in (S63l we have defined Z^{t) = 1, which we find to be a consistent analytical continuation of the GGE 
partition sum to n = 0. 

After the insertion of the string ansatz, the charges take the form Ap = Zyfi Ap{kj,mj) with 


Ai(fc,m) = km , 

A 2 {k, m) = mk^ + “ m^)<? , 

A3(fc,m) = (4fc^ + (1 - TO^)c^) . 

Ai{k, m) = mk"^ + ^c^k^m (l - m^) + (Sm"' - lOm^ + 7) 


Exchanging in (S63) the sum over Ug and the one over n, we can rewrite it as 

1 ~(n,). 


Jf}(u,t)= E — 


ns=0 


where the partition function with a fixed number of strings has been introduced 

{-uc)^i 

mi...m„.=i (27rc)"» 


Z^J'^\u,t)= E 


f 


dk 0pAp{kj,mj) 

u— --- 


0=1 


$(k, m) 


To proceed further we rewrite this expression introducing the determinant of a 71^ x matrix 

1 


det 




cmj) = d)(k, m) 


^ 2iki - 2ikj + c{mi + mj) / 
which follows by the Gauchy-determinant formula nni, and leads to 

dkj ) + Sp>i 0pAp{kj ,mj ) 


=1 


0 = 1 


27r 


det 


1 


2iki - 2ikj + c{mi + mj) ^ 


(565) 

(566) 

(567) 

(568) 

(569) 

(570) 

(571) 

(572) 


We now employ the Laplace expansion for the determinant 

1 \ 1 


det 


2iki - 2ikj + c{mi + mj) ^ 


2ikj - 2ikp. + c{mj + mp.) 

-Vj{2ikj-2ikp^-\-c{rnj+rrip^)) _ 

V 

-1 


D Jvt>0 ^ D Jvi>0 A 


-2ikj{vj-vp. )-cmj{vj+vp.) 


(S73) 


where in the last equality we exchanged P with P ^ taking into account that they have the same parity. Replacing 

det lC^{vi,Vj). 


( |S73[ ) inside ( |S72[ ) we obtain 

Z^^’\u,t) = f[ 

/ oo 

^ 2(^_^y^ll'ni‘g-tA2{k,m) + 'ZppiPpAp{k,m) ^-2ik{vx-V2)-cm{vi+V2) 

oo 27r 


with the Kernel 


Erom this it follows that 

Jp{u, t) = Det(l + Ho/Cf Ho) 

in terms of a Fredholm determinant, where Hg is the projector on [s,oo). 


(574) 

(575) 

(576) 
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MOMENTS FROM GENERATING FUNCTION 


As discussed in the text we expand the expressions for for m = 1,2,3 in the conserved charges of the model 


Ap. We display here their full expression for m = 1,2 and only the leading order in 1/n for m = 3: 


An,l — 
An ,2 = 


w( 


tiAo ~ A-\ + 




= —+ K-^ 2 ) + 0 (n) 


n{n 1 '' ' " ^ ^ 

{n^ - l)n{5n - 6)c^ (A^ - 71 ^ 2 ) 0 ^ 4(n - 1 )^ 3 ^! + (jn? - 3n + 3)A| - (n- 1 ) 71^4 + A\ - 277^2^^ 


720 


6 


(n - 3)(71 - 2) (71 - 1)77 


(S77a) 
= (S77b) 


1 1 /l2 

= — (4A1A3-A^-SA^) , + —(12A2A2 + 20AiA3-llAt-15A2-6A4)- ^c2 + 0(77) (S77c) 

677 36 6 


A, 


n,3 


12077 


[-10 {Af - 4 A 3 A 1 + 3A^) + A^ - 20 A 3 A? + 15 {Al + 2 A 4 ) A^h 


+ ( 24 A 5 - 6 OA 2 A 3 ) Al + 10 (SA^ - 6 A 4 A 2 + AAl)] + 0(77°) . (S77d) 


Using the GGE partition function in (S64), we can rewrite 

= Km{{dp})[Z^{t)] 


(S78) 


where we replaced in A„,m the charges Ap dp computed setting all /3’s to zero afterwards. We immediately see that 
the first expression (S77a) for A„,i leads to equation (|l^ in the text, valid for any 77, using the equivalence 82 = -dt 


and replacing Af which comes from the STS (see be low). Note also that in the limit 77 0 we will be allowed to 


discard in the Taylor expansion as written in ( S77a|^7d ), the terms formally of order O(n^) and higher (not written 
for 777 = 3), since (i) we have checked that they all contain derivatives acting on and (ii) Z^{t) 1 when 77 ^ 0 

(see above). 

The above expressions can be further simplified employing the STS symmetry, manifested here as the invariance 
of ( |S64[ ) under the shift of all rapidities fj,j fj,j + k for an arbitrary constant k. Such symmetry was used in e.g. 
Appendix of Ref. [T^] in presence of charges Ai, A 2 only. Here we extend the analysis to the enlarged context of the 
GGE partition sum, which contains all charges Ap. Since under the shift (Aq s 77) 


Ap + 


50 




at the order 0{k) in (S64|, we obtain the equality, valid for arbitrary /3 


-2tdi + n8i + ^p/3p9p-i 

P=2 


Z^(t) = 0. 


(S79) 


(S80) 


By expansion order by order in /3, one obtains an infinite list of equalities involving conserved charges, that can be 
generated easily by 


-2fAi + ^pAp-i^Ap 

P=i 


g({A})) = 0. 


(S81) 


where g({A}) is an arbitrary polynomial in the Ap. Here, these identities are understood as inserted in the integral (15) 


over string momenta, which is denoted here as (...). Equivalently, the charges are afterwards replaced by derivatives 
applied to Z^{t) as in (S78). Varying the choice of g({A}), we obtain for instance 


(Af^i) = 0 , 




2 kfk 


2 t 


(S82) 


where the last equality goes beyond the standard use of STS (in e.g. Ref. [E]). 

Thanks to these identities, the limit 77 -> 0 can be taken in (S77). For 777 = 1,2 only A 2 survives, that, being the 


energy, can be replaced according to 82 = - 8 t. Then using (S82), Eqs. (S77a) and (S77cl lead respectively to (16 1 and 


( 22 ) in the text, using that 


= H-dtVZnit) ^ (-5t)nnZp(0;0|f) . 


(S83) 
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For m > 3, higher conserved charges are involved and the moments are not simply related to the free-energy distribu¬ 
tion. For instance for m = 3 after using STS we find that we can replace: 


A, 


1 


n,3 


c^Al 3 Al ^Al ^Al ^Ai A2A4 


2 t 


At 


2 t 


with no further simplification. Using (S78) the final result for is thus a linear combination of 

z^{t) -1 


n^O 


(S84) 


(S85) 


/3=0 


where di-^ = dp._^. The can then be expressed as derivatives acting on the Fredholm determinant expression! S761. 

To see this, we write it as a Mellin transform 


which can be inserted in (S631 to give 


^00 ^ 

= dyZf3{y,t)y'^, 

00 ^ 

Jl3{u,t) = j^ dyZp{y,t)e~'''^ 


(S86) 


(S87) 


with Jp{0,t) = 1- This shows that Zi 3 {y,t) can be equivalently defined as the inverse Laplace transform of Jp{u,t)- 
Inserting (S86) in (|S85|) and taking the limit n 0, we obtain for k> 3 


f dyZf 3 {y,t)'^ny = .. .di^ f — {e - Jp{u,t)) = - f — (9q .. J/ 3 (u,t) 

Jo Jo u Jo u 

y 

—a,, ...9,,Det(l + no/C(^no). (S88) 

0 u 


du 
/o u 

Derivatives of a FD are easily evaluated, for instance for A: = 2 one obtains 


dad^MuA) = [Tr[(l + + Tr[(l + Tr[(l + d^lC^Y 

-Tr[(l + /C(^)-iaa/C(^(l + /Cf)"iab/Cf]]Det(l + no/C(^no) . (S89) 

The derivatives with respect to the Lagrange multiplier /3, once applied to the kernel can be converted, using 
(S75) and (S65), into a combination of derivatives with respect to vi,V 2 , e.g. 


dil^uin,V2) = 


2 ic 


8c 


o ^ ^ C^di/C^(vi,V2) 

a3/L^(l’l,'C2j = -TT- + 


64c 


and more generally, from (|S75[) 


djlC^(vi,V2) = Ap{k = '-d, 


p\^ —^Vi+V 2 ') * 


(590) 

(591) 

(592) 


These formulas can be iterated in case of multiple derivatives, leading to explicit but lengthy expressions. The limit 
/3 -> 0 can now be taken safely and we then have to deal only with the standard kernel /C„ s as given in Ref. 

[7]. The calculation at arbitrary time by this method is possible but very demanding. 

Let us now show how the asymptotics at large time simplifies. We set = YtjA and introduce s by = (uc)e~^. 
Moreover, we rescale k kc/X and Ani. 2 /c, so that lCu{vi,V 2 ) ^ XlCs{vi,V 2 )jc, because of the integration 

measure, and we obtain finally the A ^ oo limit 


/c: 


die r 

— / dyA\{y + k'^ + s + Vx + V2)e~'‘’"'''’^~'’'^'^=-2^I^Kp,^{2^l^{vx + sl2),2^l^{v2 + sl2)). (S93) 

oo 27r Jo 


with the Airy Kernel defined as 


oo 

KAi{x,y)= / dy ki{vi + y) k\{v 2 + y) 

Jo 


Ai(ui) Ai'(u 2 ) - Ai'(ni) Ai(z; 2 ) 


fl - V 2 


(S94) 
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In this regime, the Fredholm determinant can be computed efficiently numerically following m- 

Let us show how to obtain the leading large time behavior on the example of m = 3. First one sees that one can 
decompose the charges (S65) as: 


A2{k,m) = A2{k,m) + —m , 
c 2 ~ 

A3(k,m) = A3{k,m) + —Ai{k,m) , 

7c'^ - 

Ai{k,m) = A4{k,m) + —A2{k,m) 


(595) 

(596) 

(597) 


and Ai{k,m) = Ai{k,rri), where Ap(k,m) is the homogeneous part, i.e. Ap{^,m) = A Ap{kc, Xm) under the 
rescaling k kc/X. Thus in (S84) one can formally decompose the charges as: 


A 2 - A 2 + 


c^n 


A^i - A^ + 


c^Ai 


Aa — Aa + 


7c'^n c^A2 


(S98) 


12 ’ 4 ’ 240 2 

and the above analysis of the FD at large time shows that the hat charges scale as Ap ~ which we denote 

by ’’normal” scaling, and Z„/n scales as . For example A 2 Znln — —i -9tln^r;(0;0|t) = c^/12 + 0(t has 


“anomalous” scaling, while A 2 Znln °> 0{t has “normal” scaling with time at large time. 
Putting all together, in the limit n ->■ 0 we find: 




2^c®/3 


9(»/» 

where the sub-leading term is obtained using (S98l in (S84). Then the term scaling as gives at large times 


(S99) 


y 4t 


3 A 2 A^ A 




n 12t [ t 2 ) ^ ^ 


9t5/3 


(SlOO) 


CONJECTURE FOR lnp„(t) 

In [13] relations between some discrete DP models and the eigenvalues of the GUE (and LUE) random matrix 
ensembles are discussed. For instance, calling the maximum energy of an ensemble of k non-crossing paths of ^ N 
steps with nearest neighbor endpoints in a semi-discrete DP model, and Ai > A 2 the two largest eigenvalues of the 
GUE(N), it is stated that: 

hi=dXi (SlOl) 

^2 “ hi =d X 2 (S102) 

where =d denotes equality in law (where, for our purpose, we consider N large). Under appropriate rescaling |8], 
universality then suggests that at large time for the continuum DP model: 

lnZ,:2-^+^(c2i)l/3 (S103) 

(S104) 

where the limit of small e is implicit, and we denote = Zp{0;Q\t). Here X 2 > X 2 ^^^e two largest eigenvalues of 
the GUE, scaled so that the largest one obeys the Tracy Widom cumulative distribution F 2 . These equations imply: 

\npp{t) = In -2\nZr,c. « -im43{chfl^ (S105) 

with X 2 = -I.77I0868 and X 2 = -3.6754 (value given in (HI). We have obtained preliminary numerical support for this 
behavior m This means that in typical disorder configurations Pr]{t) decreases (sub)-exponentially fast with time, 
while its moments are dominated by rare, but not so rare, disorder configurations. This is similar to the behavior of 
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the probability qrj{t) that a single DP does not cross a hard wall at a; = 0 studied in [T5], where it was found that 
Ing^(t) = -(l^-~ -1.49134(c^t)^'^^ where Xi is distributed according to the GSE TW distribution with 
^ = -3.2624279. Note that avoiding a second polymer (which itself wanders and competes for the same favorable 
regions of the random potential) is more restrictive in phase space for the DP than a hard wall, consistent with 
Inp^(t) < ln( 7 ,,(t) as we find. 
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